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Thermal fluctuations and melting transitions for rotating single-component superfluids have been intensively 
studied and are well understood. In contrast, the thermal effects on vortex states for two-component superfluids 
with density-density interaction, which have a much richer variety of vortex ground states, have been much 
less studied. Here, we investigate the thermal effects on vortex matter in superfluids with U(l) x U(l) broken 
symmetries and intercomponent density-density interactions, as well as the case with a larger SU(2) broken 
symmetry obtainable from the U(l) x U(l)-symmetric case by tuning scattering lengths. In the former case 
we And that, in addition to first-order melting transitions, the system exhibits thermally driven phase transitions 
between square and hexagonal lattices. Our main result, however, concerns the case where the condensate 
exhibits SU(2)-symmetry, and where vortices are not topological. At finite temperature, the system exhibits 
effects which do not have a counter-part in single component systems. Namely, it has a state where thermally 
averaged quantities show no regular vortex lattice, yet the system retains superfluid coherence along the axis 
of rotation. In such a state, the thermal fluctuations result in transitions between different (nearly)-degenerate 
vortex states without undergoing a melting transition. Our results apply to multi-component Bose-Einstein 
condensates, and we suggest how to experimentally detect some of these unusual effects in such systems. 

PACS numbers: 67.25.dk, 67.60.Bc,67.85.Fg,67.85.Jk 


I. INTRODUCTION 

Bose-Einstein condensates (BECs) with a multicomponent 
order parameter, and the topological defects such systems sup¬ 
port, represent a topic of great current interest in condensed 
matter physics. H21 Such multicomponent condensates may be 
realized as mixtures of different atoms, mixtures of different 
isotopes of an atom, or mixtures of different hyperfine spin 
states of an atom. The interest in such condensates from a 
fundamental physics point of view is mainly attributed to the 
fact that one may tune various interaction parameters over a 
wide range in a BEC. This enables the study of a variety of 
physical effects which are not easily observed in other super¬ 
fluid systems such as He^ and He^. 

The behavior of a single-component BEC under rotation is 
well known. The ground state is a hexagonal lattice of vortex 
defects which melts to a vortex liquid via a flrst-order phase 
transition. This is well described by the London model, where 
amplitude fluctuations may be ignored. Over the years, in the 
context of studying vortex lattice melting in high-Tc super¬ 
conductors, many works have confirmed this through numeri¬ 
cal Monte Carlo simulations for systems in the frozen g auge, 
three-dimensional (3D) XY and Villain approximations,®^ 
as well as in the lowest-Landau-level approximation,!^ and 
by mapping it to a model of 2D bosons.!^ Single compo¬ 
nent cond ensat es have been available experimentally for quite 
some tim^^l^, and the hexagonal lattice ground state has been 
verified.!^ 

Condensates with two components of the order parameter 
have also been studied extensively. Analytical works focus¬ 
ing on determining the T = 0 ground states have demon¬ 
strated interesti ng vortex solu tions and a range of unusual 
lattice structureBy varying the ratio between 
inter- and intra-component couplings, the ground-state lat¬ 


tice undergoes a structural change from hexagonal symmetry 
through square symmetry to double-core lattices and interwo¬ 
ven sheets of vortices. Similar systems with three compo¬ 
nents have also been studied.C^ Experimentally, spinor con¬ 
densates have been realized in two general classes of systems. 
The first option is to use one species of atoms, usually rubid¬ 
ium, and prepare it in two separate hyperfine spin states.^ 
VorticeP and vortex latticeP have been realized in these bi¬ 
nary mixtures, where both hexagonal and square vortex lattice 
states were observed. The other o ption is to mix condensates 
of two di fferen t species of atoms.The use of Eeshbach 
resonanceallows direct tuning of the scattering lengths, 
and by extension the inter- an d intr acomponent interactions of 

multicomponent condensates. 

In this paper, we consider a specific model of a two- 
component BEC, which has the full range of fluctuations of 
the order parameter field included, as well as intercompo¬ 
nent density-density interactions. We consider the model with 
U(l) X U(l)- and as SU(2) symmetries. Eor the U(l) x U(l) 
case, we And a succession of square and hexagonal vor¬ 
tex ground-state patterns as the intercomponent interaction 
strength is varied, along with the possibility of thermal re¬ 
construction from a square to a hexagonal vortex lattice as 
temperature is reduced. 

The SU(2)-symmetric case is interesting and experimen¬ 
tally realizable. In this case U(l)-vortices are no longer topo¬ 
logical, in contrast to the U(l) x U(l)-symmetric case. In this 
case, when fluctuation effects are included we And a highly 
unusual vortex state where there is no sign of any vortex lat¬ 
tice. Nonetheless, global phase coherence persists. This state 
of vortex matter is a direct consequence of massless amplitude 
fluctuations in the order parameter, when the broken symme¬ 
try of the system is SU(2). At the SU(2) point, but at lower 
temperatures, we also observe dimerized vortex ground-state 
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patterns. 

The paper is organized as follows. The model and defini¬ 
tions of relevant quantities are presented in Section |I^ The 
technical details of the Monte Carlo simulations are briefiy 
considered in Section In Section |IV| the results are pre¬ 
sented and discussed. In Section |Vj we discuss how to ex¬ 
perimentally verify the results we find. Some technical de¬ 
tails, and the investigation of the order of the melting transi¬ 
tions with full amplitude distributions included, for the cases 

= 1 and = 2, are relegated to Appendixes. 


where rriij = rrii mj/{mi + mj) is the reduced mass. In 
this paper we focus on using BECs of homonuclear gases 
with several components in different hyperfine states, hence 
rrii = m^i. Inter-component drag in BEC mixtures has been 
considered in previous works using Monte-Carlo simulation 
(ignoring ampfitude fiuctuations), but we will not consider this 
case here.™ 

We find it convenient for our purposes to rewrite Q on 
the following form, the details of which are relegated to Ap¬ 
pendix]^ 


II. MODEL AND DEFINITIONS 

In this section we present the model used in the paper, first 
in a continuum description and then on a three-dimensional 
cubic lattice appropriate for Monte Carlo simulations. The 
relevant quantities for the discussion are also defined. 


A. Continuum model 

We consider a general Ginzburg-Landau(GL) model of an 
A^-component Bose-Einstein condensate, coupled to a uni¬ 
form external field, which in the thermodynamical limit is de¬ 
fined as 
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is the Hamiltonian. Here, the field A'^ formally appears as a 
non-fiuctuating gauge-field and parametrizes the angular ve¬ 
locity of the system. The fields are dimensionful complex 
fields, i and j are indices running from 1 to A^ denoting the 
component of the order parameter (a “color” index), a'- and 
g'^j are Ginzburg-Landau parameters, = /?^/2e is the the 
coupling constant to the rotation induced vector potential, and 
rrii is the particle mass of species i. Eor mixtures consisting of 
different atoms or different isotopes of one atom, the masses 
will depend on the index i, while for mixtures consisting of 
atoms in different hyperfine spin states, the masses are inde¬ 
pendent of i. The inter- and intracomponent coupling parame¬ 
ters g[j are related to real inter- and intracomponent scattering 
lengths, aij, in the following way 


4:7rJi?aii 

(3) 


rrii 
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H = 





(5) 


Here, is an A/'-component spinor of dimensionless com¬ 
plex fields, which consists of an amplitude and a phase, ijji = 
\^jji \ exp {iOi), — i^A'^ is the covariant derivative, 

and summation over repeated spatial indices is implied. We 
neglect, for simplicity, the presence of a trap and the centrifu¬ 
gal part of the potential. We consider only the case where the 
vector potential is applied to each component of as follows 
from the fact that the masses are independent of species-index 


We have studied this model in detail with N = 2, where we 
write the potential in the form 

V{^) -if ( 6 ) 

This formulation is more relevant for our discussion, as it im¬ 
mediately highlights the symmetry of as well as the soft 
constraints applied to it. The details of the reparametrization 
are shown in Appendix [A| 

Note that Eq. [^may also be rewritten in the form (correct 
up to an additive constant term) 

V= (7? + w)(|tAi|'‘ + |V;2|^)+2(?7-a;)|t/)i|^|?/)2|^. (7) 

Comparing with Eq. we have gn = ^22 = g = g uj and 
gi 2 = g — The model features repulsive inter-component 
interactions provided g — uo > 0, and this is the case we will 
mainly focus on. We will however briefiy touch upon the case 
g — uj < 0 corresponding to an attractive inter-component 
density-density interaction, which leads to ground states with 
overlapping vortices in components 1 and 2. Normalizability 
of the individual order parameter components, or equivalently 
boundedness from below of the free energy, requires that g A 
cc > 0. Thus, while uj > g makes physical sense, uo < —g 
does not. In this paper, we assume 7 ^ > 0 and cc > 0. 

Two-component BECs feature considerably richer physics 
than a single-component BEC. Since the gauge-field 
parametrizing the rotation of the system is non-fiuctuating, 
there is no gauge-field-induced current-current interaction be¬ 
tween the two condensates (unlike in multi-component su¬ 
perconductors). The only manner in which the two super- 
fluid condensates interact is via the inter-component density- 
density interaction 2{g — uj)\2pi\‘^\2p2\‘^ - In the limit where 
the amplitudes of each individual component are completely 
frozen and uniform throughout the system, one recovers the 
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physics of two decoupled 3DXY models, with a global 
U(l) X U(l) symmetry. The density-density interaction be¬ 
tween jjji and ^2 leads to interactions between the topological 
defects excited in each component. As a result, a first order 
melting of two decoupled hexagonal lattices is not the only 
possible phenomenon that could take place. Previous experi¬ 
ments and numerical studies have reported a structural change 
of the ground state from a hexagonal to a square lattice of vor¬ 
tices as the effective inter-component coupling is increased.!^ 
This corresponds to increasing the ratio r]/uj m our case. As 
we will see below, other unusual phenomena can also occur, 
notably when thermal fiuctuations are included. 

One special case of the model deserves some extra atten¬ 
tion. If one takes the limit cc ^ 0 in Eq. © the symmetry 
of the model is expanded to a global SU(2) symmetry. One 
may then shift densities from one component to the other with 
impunity, as long as is left unchanged. This ef¬ 

fectively leads to massless amplitude-fluctuations in the com¬ 
ponents of the order parameter. Therefore, it is possible to 
unwind a 27r phase winding in one component by letting the 
amplitude of the same component vanish. The introduction 
of this higher symmetry leads to very different vortex ground 
states than what are found in the U(l) x U(l)-symmetric case 
with cc 7 ^ 0. 


B. Separation of variables 

In multicomponent GL models complex objects, such as 
combinations of vortices of different colors, are often of in¬ 
terest. In general, it is possible to rewrite an A^-component 
model coupled to a gauge field, fiuctuating or not, in ter ms of 
one mode coupled to the field and — 1 neutral modes. 

For a more general discussion of charged and neutral modes 
in the presence of amplitude fiuctuations see Refs. [42] and 
l43l Considering only the kinetic part of the two-component 
Hamiltionian, we have the following expression: 




2|T^|' 
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Hence, the first mode couples to the applied rotation, while the 
second does not. This corresponds to the phase combinations 
0i -f- O 2 and 0i — O 2 , respectively. 


C. Lattice regularization 


and introducing real phases and amplitudes \ " 

we can rewrite the Hamiltonian. 


H ^ ^ |'0r+/x,i I I COs(^r-|-^^2 ^/x,r) 

r,A 

i 

+ J2Vi^r). ( 10 ) 

r 

The lattice spacing is chosen so that it is smaller than the 
relevant length scale of variations of the amplitudes. A di¬ 
mensionless vector potential, has also been introduced. 
See Appendix [A] for details. We denote the argument of the 
cosine as x^-, as a shorthand. 


D. Observables 


An important and accessible quantity when exploring phase 
transitions is the specific heat of the system. 


cv = 


(H^) - {Hf 
£3 


( 11 ) 


While crossing a first order transition there is some amount 
of latent heat in the system, manifesting itself as a (^-function 
peak of the specific heat in the thermodynamic limit. On the 
lattice one expects to see a sharp peak, or anomaly, at the tran¬ 
sition. This is used to characterize the transition as first order. 

A useful measure of the global phase coherence of the sys¬ 
tem, is the helicity modulus, which is proportional to the su- 
perfiuid density. It serves as a probe of the transition from 
a superfiuid to a normal fluid. In the disordered phase, the 
moduli in all directions are zero, characterizing an isotropic 
normal-fiuid phase. The cause of this is a vortex loop blowout. 
Moving to the ordered phase, all moduli evolve to a finite 
value. If we turn on the external field we still have zero coher¬ 
ence in all directions in the disordered phase. In the ordered 
phase, however, the helicity modulus along the direction of 
the applied rotation jumps to the finite value through a first or¬ 
der transition. The value of the transverse moduli will remain 
zero. Formally, the helicity modulus is defined as a deriva¬ 
tive of the free energy with respect to a general, infinitesimal 
phase twist along r^,®. That is, we perform the replacement 


In order to perform simulations of the continuum model, 
we define the field on a discrete set of coordinates, i.e 
^(r) ^ where r G (ix -j- jy -1- /cz|i, j, /c = 

Here, L is the linear size in all dimensions; the system size 
is V = L^. We use periodic boundary conditions in all di¬ 
rections. By replacing the differential operator by a gauge- 
invariant forward difference 


^r,i —^ i — 

in the free energy, and calculate 


T 




d^F[0'] 

d^l 



( 12 ) 


(13) 
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Here, b = ( 6 i, 62 ) represents some combination of the phases 
0i and O 2 , biOi + 626 ^ 2 . To probe the individual moduli, bi is 
chosen as bi = (1,0) or bi = (0,1). The composite, phase 
sum variable is represented by the choice bi = ( 1 , 1 ), while 
bi = (1, —1) is the phase difference. Generally, for a two- 
component model, the helicity modulus can be writte n as the 
sum of two indivudual moduli, and a cross term,*^^^ 


+ ^2^Ai,(0,l) + 26 i62T^,12. (14) 

For the model considered in this paper, the individual helicity 
moduli can be written as 




X]V’r,V’r+A.i COs(xif,i) 

r 

- /?( (y]^r,iV’r+A.i sin(Xr,i) 


while the mixed term has the form 


(T^,12) = -^/(y^lV’r.il IV’r+A,i|sin(xJf,i) 

\ \ r 

y]|V’r,2| |V’r+A,2|sin(xJf,2) 


(15) 


(16) 


We denote the helicity modulus of the phase sum (i i) as 
T+ as a short hand. 

fj, 

The structure factor, 5'i(q±), can be used to determine the 
underlying symmetry of the vortex lattice. Square and hexag¬ 
onal vortex structures will manifest themselves as four or six 
sharp Bragg peaks in reciprocal space. In a vortex liquid phase 
one expects a completely isotropic structure factor. The struc¬ 
ture factor is defined as the Fourier transform of the longi¬ 
tudinally averaged vortex density, {ni{T±)), which is subse¬ 
quently thermally averaged. 





(17) 


Here (r^) is the denisity of vortices of color i averaged over 
the ^-direction 


ni{r±) = (18) 

z 

and r± is r projected onto a layer of the system with a given 
^(-coordinate. The vortex density is calculated by traversing 
each plaquette of the lattice, adding the factor Xir 
link. Each time we have to add (or subtract) a factor of 27r 
in order to bring this sum back into the primary interval of 
(—TT, tt] a vortex of color i and charge -Hl(-l) is added to this 
plaquette. 

In addition to the structure factor, we look at thermally av¬ 
eraged vortex densities, {ni{r±)), as well as thermally and 
longitudinally averaged amplitude densities, {\i^if (i*±)), de¬ 
fined similarly to Eq. ( p^ . 


|V’i|^(r±) = (19) 

z 

This provides an overview of the real space configuration of 
the system. 

When including amplitude fiuctuations, which, when the 
potential term is disregarded, are unbounded from above, it 
is of great importance to make sure all energetically allowed 
configurations are included. To this end, we measured the 
probability distribution of I'lpif, P{\^pif) during the simula¬ 
tions by making a histogram of all field configurations at each 
measure step, and normalizing its underlying area to unity in 
post-processing. 

The uniform rotation applied to the condensates is imple¬ 
mented in the Landau gauge: 


A = (0,27r/x,0), (20) 

where / is the density of vortices in a single layer. Note that 
this implies a constraint Lf G (1,2,3,...) due to the periodic 
boundary conditions. When probing a first order melting tran¬ 
sition, it is important to choose a filling fraction large enough 
that an anomaly in the specific heat is detectable. However, 
if the filling fraction is too large, one may transition directly 
from a vortex liquid into a pinned solid, completely missing 
the floating solid phase of interest. This scenario is character¬ 
ized by a sharp jump in not only the longitudinal, but also the 
transverse helicity modulus .1^^^ One must therefore chose / 
small enough, to assure that the vortex line lattice is in a fioat- 
ing solid phase when it melts. 


III. DETAILS OF THE MONTE CARLO SIMULATIONS 

The simulations were performed using the Metropolis- 
Hustings algorithm.1^3111 phase angles were defined as G 
(— 7 r, 7 r], and amplitudes as G (0,1 + dtl)]. The choice 
of (5^/) will be discussed further, as it is important to ensure 
inclusion of the full spectrum of fiuctuations. Both the phases 
and the amplitudes were discretized to allow the use of ta¬ 
bles for trigonometric and square root functions in order to 
speed up computations. We typically simulated systems of 
size = 64^, with sizes up to = 128^ used to resolve 
anomalies in the specific heat. We used 10^ Monte Carlo 
sweeps per inverse temperature step, and up to 10 ^ close to the 
transition. 10 ^ additional sweeps were typically used to ther- 
malize the system. In the simulations, we examined time se¬ 
ries of the internal energies taken during both the thermaliza- 
tion runs and the measurements runs to make sure the simula¬ 
tion converged. One sweep consists of picking a new random 
configuration for each of the four field variables separately in 
succession, at each lattice site. Measurements were usually 
performed with a period of 100 sweeps, in order to avoid cor¬ 
relations. Ferrenberg-Swendsen multi-histogram reweighting 
was used to improve statistics around simulated data points, 
and jackknife estimates of the errors are used. 
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FIG. 1: (Color online) The probability distribution of the 
amplitudes, for = 2 , at inverse temperature 

P = 1.20, / = 1/32, and f] = 2, with uj values from 0 to 3. 
The distribution is completely symmetric in i. 


Fig. p shows the probability distribution of the amplitudes. 
We get a peaked distribution for finite cj. On the 
other hand, when cc = 0, this is no longer the case. The dis¬ 
tribution now approaches a uniform distribution on the inter¬ 
val (0,1]. In this case the parameter r] serves to control the 
approach to uniformity, r] ^ oc corresponding to the CP^ 
limit. 

With these initial simulation runs as a basis, we choose 6'ip 
appropriately in order to capture the entire spectrum of fluctu¬ 
ations. 


IV. RESULTS OF THE MONTE CARLO SIMULATIONS 


In this section, the — u; phase diagram of ground states 
is explored by slow cooling and examination of vortex and 
amplitude densities, as well as structure factors. In addition 
to the expected hexagonal and square vortex ground states, 
several interesting regions of the parameter space are inves¬ 
tigated further. A special case between the square and the 
hexagonal region of the phase diagram is discovered, where 
the lattice first forms a square structure, but thermally recon¬ 
structs into a hexagonal lattice as the temperature is decreased 
further. Furthermore, we consider in detail the u; = 0 line in 
the phase diagram, where we discover additional vortex fluc¬ 
tuation effects. For cc = 0, the system features an SU(2) sym¬ 
metry. An unusual feature is an interesting state with global 
phase coherence, but without a regular vortex lattice. In this 
case ordinary vortices do not have topological character due to 
SU( 2 ) symmetry. Additionally, we obtain several interesting 
vortex structures characterized by dimer-like configurations at 
lower temperatures. Here, we observe honeycomb lattices, or 
double-core lattices, and stripe configurations, consistent with 
previous T = 0 results.^ 

We also examine the melting transitions of the square and 


hexagonal lattices with the full amplitude distribution in¬ 
cluded, as well as the melting of the hexagonal lattice in a 
model with = 1 as a benchmark of the method. To classify 
the transition, we look at thermal averages of the specific heat, 
helicity moduli, and vortex structure factors. These results are 
presented in Appendixes [B] and [C| 


A. The r] — (jj phase diagram 


Adding a second matter field and inter-component density- 
density interactions results in a considerably richer set of 
ground states than in the single-component case. In the ab¬ 
sence of a fluctuating part of the rotational ’’gauge-field” there 
will be no gauge-field-mediated inter-component current- 
current interactions. For 77 — cc < 0, (? 7 , cc) >0) the effective 
inter-component density-density coupling r] — uj is negative 
and the ground state of each color of condensate has a hexag¬ 
onal symmetry, as shown in Fig. |2(a)[ If, on the other hand 
T] — UJ > 0 , the inter-component coupling becomes positive. 
Now, for sufficiently large ratios r]/uj, the vortices arrange 
themsel ves in to two inter-penetrating square lattices, shown 
in Fig. 2(b) The value of the ratio r]/uj for which the lat¬ 
tice reconstructs depends on the strength of the rotation, /. If 
we neglect fluctuations, 77 — cc < 0 is expected to result in a 
hexagonal lattice, while 77 — cc > 0 leads to a square lattice for 
sufficiently large r]/uj. 

The physics of the reconstruction of the lattice can be ex¬ 
plained by modulations of the amplitude fields. The existence 
of static periodic amplitude-modulations (density-variations) 
is due to the presence of vortices. Without vortices (/ = 0) 
and UJ > 0 , the ground state is one where both amplitudes are 
equal and smooth. Vortices in one component tend to suppress 
locally the corresponding amplitude, which in turn means that 
the term 77 ( 17/1 p -h I 7 / 2 P — 1 )^ enhances the amplitude of 
the other component. At small uj, i.e. large r] — uj there is 
a strong tendency to form a square density lattice due to this 
intercomponent density-density interaction. Conversely, if uj 
is large enough compared to 77 , the density-density interac¬ 
tion is not strong enough to overcome the isotropic current- 
current interactions between same-species vortices. In other 
words if the current-current interactions dominate the inter¬ 
species density-density interactions, a hexagonal lattice is en¬ 
ergetically favoured over a square lattice, and vice versa. 
Note that similarly a square vortex lattice forms in two- 
component London models with dissipationless drag when 
there are competi ng in ter- and intra-species current-current 
vortex interactions .EEH 

Figures. [3(a)| and [3(b)] show the phase diagrams for filling 
fractions / = 1/32 and / = 1/64, respectively. The sep¬ 
aration line is approximate and drawn from several separate 
simulations. 

To clarify what is going in Figs. |3(a)] and |3(b)[ we refer to 
Figs. [TT] and in Appendix Here, we show tableaus to 
illustrate in more detail how the density and vortex lattices 
reconstruct at a temperature well below any melting tempera¬ 
tures of the vortex (and density) lattices, as the density-density 
interaction 2(77 — ci;)| 7 /ip| 7 / 2 p is varied. Specifically, we fix 
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FIG. 2: (Color online) Representative configurations of the 
two main ordered phases in the U(l) x U(l) region, (a) 
shows a square structure at ( 77 , cc) = (5.0,0.5), while (b) 
illustrates the hexagonal structure at ( 77 , cc) = (5.0, 5.0). 

Each subfigure shows vortex densities, {ni{T±), in the left 
column, amplitude densities, {I'lpif (i*±), in the right column, 
and structurue factors (insets) of each component as 
indicated. The induced vortex density and inverse 
temperature are fixed to / = 1/64 and p = 1.5 in both 
subfigures. 


the interaction parameter 77, as well as the inverse tempera¬ 
ture p and filling fraction /, while increasing the parameter 
uj. This reduces the effective inter-component density-density 
interaction which favors a square lattice, until the lattice re¬ 
constructs from square to hexagonal symmetry. 

When 77 = cc, it is seen from Eq. 0 that the two compo¬ 
nents of the order parameter decouple. Eor cc < 77 the inter¬ 
component density-density interaction is repulsive, while it is 
attractive for cc > 77. Eor cc < 77, the vortex lattices (and the 
density lattices) are intercalated, while for cc > 77 they are co- 



n 



n 


EIG. 3: (Color online) The 77 — cc phase diagram of the 
ground states for / = 1/32 (top) and / = 1/64 (bottom). 
The simulations were performed for a range of ( 77 , cc) pairs to 
determine the zero temperature ground state. Approximate 
demarcation lines for the phase boundaries separating 
hexagonal lattices, square lattices, and dimerized phases, 
were drawn from these results (solid lines). I denotes the 
phase where the hexagonal vortex lattices in the two 
components are cocentric, II denotes the case where the 
hexagonal lattices are intercalated, while III denotes the 
square lattice phase. The dotted line is the line cc = 77 at 
which the intercomponent density-density interaction 
2(77 — cc)1 7/1 Pi 7 / 2 P changes sign. See also Eigs, 

Appendix [D| 
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and 
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m 


centric. In Eigs. andwe illustrate the demarcation line 
between the two situations as a dotted line in the hexagonal 
phase. 

Beyond the square and hexagonal lattices we also observe 
dimer configurations of vortices for cc = 0 , which will be 
discussed further below. The calculations are consistent with 
the ground states obtained in Refs. IblandfTOl 


B. Thermally induced reconstruction of vortex lattices 

Now we move to discussion of the effects of thermal fiuc- 
tuations in these systems. Eig. shows the vortex-densities 
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in component 1 in reciprocal space, as [3 is increased, i.e., as 
temperature is reduced, in a temperature range below where 
the lattice melts. The actual melting of the two component 
lattice is discussed in Section We fix the filling fraction 
/ = 1 /64, as well as the interaction parameters rj = 2 and 
uj = 0.5. 

For the highest temperatures shown in Fig.j^the vortex lat¬ 
tice is square. Upon cooling the system, the vortex lattice re¬ 
constructs into a hexagonal lattice, consistent with the ground 
state phase-diagram of Fig. j3(b)| The density-density interac¬ 
tion term 2 (77 — cc) | p 1 7 /^ 2 1 ^ aids formation of a square lattice 
at higher temperatures, while the current-current interactions 
drives the lattice towards a hexagonal configuration when it is 
cooled further. This means that the free energy per vortex of 
the square lattice, which is lower than that of the hexagonal 
lattice at /3 = 0.90, has become larger than that of the hexago¬ 
nal lattice when p = 1.50. This is essentially the combination 
of an energetic and an entropic effect. We observe this recon¬ 
struction not too far away from the demarcation line separat¬ 
ing a square and a hexagonal vortex lattice. Deep inside the 
hexagonal phase in Fig. |3(b)[ we observe direct vortex lattice 
melting from a hexagonal lattice to a vortex liquid. We note 
that intermediate entropically-stabilized vortex lattice phases 
were of a subject of interesting investigation in the different 
system of U{1) x U{1) superconductorhowever the vor¬ 
tex interaction form is different in this case. 


C. SU(2) vortex states 

The limit cc ^ 0 is quite different from the U(l) x U(l)- 
symmetric case uj ^ 0. From Eq. ([^, it is seen that the 
Hamiltonian is invariant under SU(2) transformations of 
Vortices, which are topological in a U(l) x U(l) model, are 
no longer topological in the SU( 2 ) case. One may unwind a 
27r phase winding by entirely transferring density of one com¬ 
ponent to the other, which may be done at zero energy cost. 

Fig. shows one of the main results of our paper. These 
are simulations with SU( 2 ) symmetry, i.e., cc = 0, as well 
as 77 = 5.0 and / = 1/64. The top panel show the phase 
stiffness associated with the phase sum, T+. This is the phys¬ 
ically relevant phase variable in this case, as it couples to the 
rotation. 

We observe that the stiffness along the ^-direction becomes 
finite at an inverse temperature, p ^ 0.9. This is what one 
would expect when a vortex lattice forms. However, the bot¬ 
tom panels, which shows the vortex density of component 1 
at /3 = 0.94, shows no apparent signs of vortex ordering. 
Hence, we have an unusual situation. There is a relatively 
large /3-range where we have a finite z-directed helicity mod¬ 
ulus of the phase sum, but no apparent ordering of induced 
vortices. A finite helicity modulus generally means that there 
are straight vortex lines with very little transverse fiuctuations 
threading the entire system along the direction in question. In 
the U(l) picture this corresponds to a regular vortex lattice. 
For an SU(2)-condensate, this is no longer the case. Large 
relative amplitude fiuctuations can occur since they have zero 
energy cost in the ground state as the energy is no longer min¬ 
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(b) 



q^ {2'n/L^a) q^ {2'n/L^a) 
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FIG. 4: (Color online) Thermally induced reconstruction 
from a square vortex lattice in either of the components at 
77 = 2,cc = 0.5, toa hexagonal vortex lattice, as /3 is 
increased. Here, / = 1/64. (a)-(f) show inverse temperatures 
P = {0.80,0.90,1.20,1.30,1.34,1.38}, respectively Each 
subfigure shows 5'i(qx) only; 5'2(q±) is identical. The 
physical reason for the reconstruction originates with the 
inter-component density-density interaction term 
2(77 — cc) It/i P 1^/2 P), and is explained in detail in the text. 


imized by a preferential value of |7/i |^ — I 7 / 2 1^- This results in 
many (nearly) degenerate vortex states between which the sys¬ 
tem can fluctuate, thus greatly simplifying the effort of mov¬ 
ing an entire, almost straight, vortex line. We are left with 
a phase where we have coherence along the ^-direction, but 
no regular vortex lattice appears in thermal averages. Nearly 
straight vortex lines will shift between a large number of de¬ 
generate, or nearly degenerate, states at a time scale shorter 
than a typical Monte-Carlo run. 

The bottom panels of Fig.|^show some inhomogeneities of 
the vortex densities, exemplifying that this is not an ordinary 
vortex liquid with segments of vortex lines executing trans¬ 
verse meanderings along their direction, which would yield 
zero helicity modulus along the direction of the field-induced 
vortices. Rather, what we have is a superposition of many 
lattice-like states of nearly straight vortex lines, where the 
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FIG. 5: (Color online) Illustration of the observed state with 
coherence along the direction of the rotation axis without a 
regular vortex lattice, seen only with SU(2) symmetry. The 
parameters used are cj = 0.0, A = 5, and / = 1/64. The top 
panel shows the helicity modulus of the phase sum, T+. The 
two bottom panels show the vortex densities, ni(r^), at 
/3 = 0.94, where the z-directed modulus is clearly finite. The 
bottom left and bottom right panels are taken from 
simulations using 10^ and 10^ Monte-Carlo sweeps, 
respectively. No apparent vortex line structure is seen here, 
and by increasing the number of Monte-Carlo sweeps the 
variations of the vortex density are smoothed out further. 
Note how the value of the average vortex density seems to 
converge towards 1/64. 
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FIG. 6 : (Color online) Two examples of SU(2) vortex 
configurations from a single simulation, for two different 
inverse temperatures. The parameters r/, /, and cj are fixed in 
each subfigure, at 77 = 1.0, cj = 0.0, and / = 1/64. (a) 
shows p = 0.84, while (b) shows p = 1.50. Each subfigure 
shows vortex densities, {ni{Y_\_), in the left column, 
amplitude densities, the right column, and 

structure factors (insets) of each component as indicated. 
This illustrates the degeneracy of the vortex line lattice in the 
isotropic limit, as the configurations evolve when P is varied. 

See appendix l^for more details. 


fiuctuations are largely collective excitations of entire nearly 
straight lines, rather than fiuctuations of smaller segments of 
lines. 

We emphasize again that these collective excitations orig¬ 
inate with large amplitude fiuctuations due to the SU{2)- 
softness of the amplitudes of the components of the superfiuid 
order parameter, rather than with phase fiuctuations. Increas¬ 
ing the number of Monte-Carlo sweeps by an order of magni¬ 


tude smooths these variations out (without noticably alteri^ 
the value of T+), as seen in the bottom right panel of Fig. 
Note how the average value of the vortex density seems to 
converge towards 1/64. This is what we expect for a vortex 
lattice or liquid in a U(l) x U(l) symmetric model, as the 
density of thermal vortices will average to zero, and / is the 
average flux density per plaquette. 

As the system is cooled further, the movements of large 
vortex lines cease, and a regular vortex lattice appears. How- 
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ever, degeneracy must still be present, as the exact pattern 
formed by the lattice is distinctively different between sim¬ 
ulations (keeping all parameters equal). The lattice also has 
a tendency to shift between configurations as the temperature 
is varied, below the temperature of initial vortex lattice for¬ 
mation. We observe two distinct classes of vortex states, illus¬ 
trated in Fig. [6| Th e two are stripes (Fig. |6(a)| ) and honeycomb 
lattices (Fig. [6(b)| ), both of which are seen in Ref.[6l Note that 
these vortex densities are taken from a single simulation, after 
the lattice has formed. Within the accuracy of our simulations 
the obtained states are not metastable. The evidence of this 
is obtained by performing several independent runs from dif¬ 
ferent initial configurations. Again, we refer to Appendix [D| 
where Fig. [^illustrates the degeneracy in the vortex line lat¬ 
tices obtained in the isotropic limit in further detail. 


V. EXPERIMENTAL CONSIDERATIONS 


Hexagonal and square lattices have already been observed 
in binary condensates of rubidium.l^ However, an SU(2) con¬ 
densate has not been realized experimentally. In this section, 
we briefiy outline under what circumstances an observation of 
an SU(2) vortex state may be feasible. 

In order to experimentally realize SU(2) conditions, one 
requires a two-component EEC, where both intra- and inter¬ 
component interactions are equal. As we have seen, the SU(2) 
physics crucially depends on this, since even minor deviations 
from this condition immediately yield U{1) x U{1) physics. 
This corresponds to cc = 0 in our parametrization. Intra- 
and inter-component density-density interactions are given in 
terms of scattering lengths. While tuning of these in an exper¬ 
iment is possible with Feshbach resonances, it may still be a 
challenge to tune two scattering lengths independently to be 
equal to a third, to arrive at the SU(2) point. From what is 
known for scattering lengths of real systems, it appears that a 
mixture of two species of the same atom, but in different hy- 
perfine states, lends itself more readily to a realization of an 
SU(2) condensate than a mixture of different atoms or a mix¬ 
ture of different isotopes of the same atom. This is so, since 
in the former case, the relevant scattering lengths typically a 
priori are much more similar to each other than what they are 
in more heterogeneous mixtures. 

One promising candidate therefore appears to be a con¬ 
densate of ^^Rb prepared in the two hyperfine states \F = 
l,m/ = 1) = |1) and \F = 2,m/ = —1) = |2). In this 
system, the three relevant s-wave scattering lengths already 
have values close to the point of interest, an = 100.da^, 
022 = 95. 00a and ai 2 = OT.fifia^, where a^ is the Bohr 
radius.Reference fTT] reports on a magnetic Feshbach res¬ 
onance at a field of approximately 9.1 Gauss, where control 
of ai 2 of the order of lOa^ is possible. Additionally, Ref¬ 
erence [51] reports on an optical Feshbach resonance of the 
state \F = l^rriF = —1), able to tune the intracomponent 
scattering length, using two Raman lasers, with detuning pa¬ 
rameters approximately given by Ai = x 75MHz and 
A2 = 27 r X 20MHz. Here, varying A 2 tunes the value of the 
scattering length around the Feshbach resonance, while vary¬ 


ing Ai changes the width. Hence, greater control of the res¬ 
onance is possible with an optical Feshbach-resonance com¬ 
pared to a magnetic one. Presumably, there should exist opti¬ 
cal Feshbach resonances able to tune the scattering length of 
either the |1) or the |2) state, for instance the one reported to 
exist at 1007G for the 11) state .E2 This resonance should be far 
enough away from the inter-component resonance at 9.1 G to 
not cause any interference. 

This suggests one possible setup. Namely, prepare a two- 
component condensate of ^^Rb in the |1) and |2) states under 
rotation, and tune ai 2 to 022 using a magnetic field. Then, 
tune ail to the same value using optical techniques, while 
taking time-of-fiight images of the condensate. The predic¬ 
tion is that as the system is tuned through the optical Feshbach 
resonance, one should observe a hexagonal composite vortex 
lattice at sub-resonance frequencies, the non unique vortex or¬ 
dering pattern, discussed above, at a frequency where all scat¬ 
tering lengths are equal, close to the optical Feshbach reso¬ 
nance, and finally the reappearance of a hexagonal vortex lat¬ 
tice at frequencies above the frequency where all lengths are 
equal Fig.|^ The observation of a featureless rotating conden¬ 
sate would be a direct manifestation of the loss of topological 
character of f/(l)-vortices in the SU(2)-symmetric case. It 
would be interesting to study the dynamics of the vortex lat¬ 
tice in this case with methods like those used irP^. For other 
discussions of S\J(N) models in cold atoms see Refs. [54] and 
[55] 

In actual experiments, a magnetic trap is used to confine 
the condensate in a given lateral region. The effect of this on 
thermal fluctuations in vortex matter has been studied in de¬ 
tail in previous theoretical wor ks fo r the one-component case, 
without amplitude fiuctuationJ^^Ell. The effect of the trap is 
to yield a maximum overall condensate density at the center 
of the trap, while depleting it towards the edge of the trap. 
As a result, the lattice melts more easily near the edge of the 
trap. As can be inferred from the work on single-component 
meltin^^^^, the results of the present paper, where no inho¬ 
mogeneity due to a magnetic trap has been accounted for, is 
therefore most relevant to the region close to the center of the 
trap. 


VI. CONCLUSIONS 


In this paper, we have investigated a two-component 
U(l) X U(l) and SU(2) Bose-Einstein condensate with 
density-density interaction under rotation at finite tempera¬ 
ture, thereby extending previous works which calculated the 
zero-temperature ground state numerically. In the U(l) x U(l) 
case we report that thermal fluctuations can lead to a phase 
transition between hexagonal and square vortex lattices with 
increased temperature. 

In the isotropic, SU(2), limit, we have observed an inter¬ 
mediate state of global phase coherence without an accompa¬ 
nying vortex lattice in the thermally averaged measurements. 
In addition, we observe a variety of dimerized vortex states, 
such as dimerized stripes and honeycomb-like lattices, which 
exist for a wide range of temperature. These lattices could 
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be observed in binary Bose-Einstein condensates in two sep¬ 
arate hyperfine states, by precisely tuning the inter- and intra¬ 
component scattering lengths to the SU(2) point through the 
use of Feshbach resonances. 
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Appendix A: Rewriting the general Hamiltionian 


Here we present the details of rewriting Eq. into Eq|^ 
which is more suited for our purposes. We repeat the starting 
point here for convenience. 
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H= d^r 
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which on the lattice reads 
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where a is the lattice constant, and we have introduced 


A Af 


(A7) 


Next, we specialize to the case N = 2, a\ = Qii = 
922 = 9, and mi = m 2 , and define to be equal to 
Ji?/mao, which sets our length scale. Note that it should not 
be confused with the coherence length in the multi-component 
case without intercomponent density-density interaction. For 
the definition of coherence lengths in the presence of multi¬ 
ple components and inter-component density-density interac¬ 
tions, see Refs. [53and[59l The energy scale is defined as Jq 
as follows 


Jo = 


9o 


(A8) 


The coupling parameters r] and u used in this paper are de¬ 
fined by comparing the potential term of Eq. ( |A6| ) to the 
form where the soft constraints \ipif y \ip 2 \‘^ = 1 and 
= 0 are implemented. Thus, we have 

yW = vilH" + IH" - + (A9) 


with 


First, we scale the field variables and Ginzburg-Landau pa¬ 
rameters, to obtain some dimensionless quantities. 


a'- = aoai, 

a'ij = gogij, 
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This gives us the Hamiltonian, 
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The lattice version of the Hamiltionian reads 
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This model will then have the following continuum form. 
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Appendix B: First order lattice melting for iV = 1 reconsidered 

As a benchmark on simulations with amplitude fluctuations 
included, we verify the well-established flrst-order melting 
transition on this model with only a single component of the 
order-parameter held, in the presence of amplitude fluctua¬ 
tions. The added feature of the computation is that the com¬ 
plete amplitude-distribution function was utilized, through the 
methods described in Section ||n| In this case, the term in the 
potential proportional to uj in Eq. [^is absent, and the potential 
reduces to 




(Bl) 


With amplitude fluctuations neglected, this model reduces to 
the much studied uniformly frustrated 3DXY model, with 
well known results as mentioned in the Introduction of the 
paper. The model features a flrst-order phase transition mani¬ 
fested as a melting of the frustration-induced hexagonal lattice 
of vortices. [i^C3E3xhe fluctuations responsible for driving this 
transition are massless transverse phase fluctuations of the or¬ 
der parameter. 

The simulations were performed with rj = 10. Fig. |7(a)| 
shows the speciflc heat, which has strong signs of an anomaly 
at = 0.751. Fig. 7(b) shows that the anomaly in the speciflc 
heat is accompanied by a relatively sharp jump in the helicity 
modulus in the ^-direction. It is also important to note that 
the helicity moduli in the transverse directions remain zero 
throughout the transition. This indicates that the vortex lattice 
melts in a genuine phase transition, and not as a result of ther¬ 
mal depinning from the underlying numerical lattice. This is 
therefore a strong indication of a flrst-order melting transition. 
Fig. shows the vortex density and structure factor immedi¬ 
ately before and after the transition. The high-temperature 
side shows an incoherent vortex liquid, characterized by a cir¬ 
cular structure factor. The-low temperature side shows that a 
clear hexagonal structure is established as soon as the liquid 
freezes. 

We emphasize that these results are not unexpected. The 
purpose of including them here, is to demonstrate that the 
method of including amplitude fluctuations into the computa¬ 
tion of the vortex lattice melting reproduces the known result 
for N = 1 , prev iously obtained in the absence of amplitude 
fluctuationiSmHUl 13 ^^ generally believed to be correct also 
when amplitude fluctuations are included. 


Appendix C: First order lattice melting for N — 2 

We next consider the melting transition for A" = 2, 
where the intercomponent density-density interaction 2{r] — 
uj) I^^ 2 1^-term in the potential energy in Eq. comes 
into play. We consider the f/(l) x U (l)-symmetric case, i.e. 
cc 7 ^ 0. Again, the full spectrum of amplitude fluctuations is 
included, using the methods described in Section [nl| 

For parameters = (0.5,1.0,1/16) and = 

(2.0,1.0,1/16) the lattices are clearly hexagonal and square, 
respectively. The hexagonal lattice obtained for 77 = 0.5 and 




/?(iAo) 


FIG. 7: (Color online) Speciflc heat (top) and helicity moduli 
(bottom) for A = 1, / = 1/16, and 77 = 10. At /3 = 0.751 
we see a clear anomaly in the speciflc heat accompanied by a 
sharp jump in the longitudinal helicity modulus. The 
transverse moduli remain zero throughout the transition. 



FIG. 8 : (Color online) Vortex density, n(r^), and structure 
factor, 5'(q±), (inset) for A = 1, / = 1/16, and 77 = 10, at 
inverse temperatures /3 = 0.749 (left) and (3 = 0.752 (right). 
This corresponds to temperatures slightly higher and lower, 
respectively, than the transition point, /3 = 0.751. 
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cj = 1.0 was found to have a melting transition at /3 ^ 0.53. 
Fig. |9(a)| shows the specific heat with a delta-function like 
anomaly at this temperature. Around this point, we have used 
a closely spaced set of temperatures, in order to get a proper 
resolution of the anomaly. Fig. |9(b) shows the helicity moduli 
of both components. Both of the 2 ;-directed stiffnesses have a 
zero expectation value in the disordered phase, indicating no 
phase coherence. In the ordered phase, both of {Tz^i) develop 
finite expectation values which means that the system has su- 
perfiuidic properties along the direction of rotation. The two 
phases are divided by a sharp jump in the longitudinal phase 
stiffness, a characteristic of a first-order transition. The drop 
is even sharper than what was obtained for the = 1 case, 
indicating an even larger latent heat associated with the tran¬ 
sition. The X- and y- directed stiffnesses remain zero in the 
ordered state, which rules out any possibility of numerical 
pinning effects. 03^ Looking further at the insets of Fig. 9a 


which show the structure factors in the disordered and ordered 
phase, we see clear evidence of an incoherent vortex liquid at 
P < 0.53 in the left inset, while the right inset shows an or¬ 
dered hexagonal vortex liquid lattice at /3 > 0.53. 

Turning to the square lattice, now the parameters in ques¬ 
tion are 7^ = 2.0 and uj = 1.0. The transition point is located 
at /3 ^ 1.11. Fig. 10(a) shows the specific heat. Again, an 
anomaly is located at the transition point. The helicity mod¬ 
uli, shown in Fig. 110(5)] also show first order behavior. Both 
z-directed components are zero on the high temperature side, 
and develop a finite value through a sharp jump at the low 
temperature side. It is important to also consider the trans¬ 
verse components. Both {Tx,i) and are zero through¬ 

out the area of interest. Here we note that the x-directed mod¬ 
ulus drops to a tiny negative value at a point after the transi¬ 
tion. This is a non-physical effect, most likely caused by a 
metastable state. We believe this is simply a numerical arti¬ 
fact, as we used a lower amount of Monte Carlo time away 
from the transition. Turning our attention to the structure fac¬ 
tors, shown in the insets of Fig. 10b we again see the isotropic 
vortex liquid in the disordered side of the transition, the or¬ 
dered side shows a square four-fold symmetry. 

Thus, both the square and the hexagonal lattices undergo 
first order melting transitions from their respective ordered 
phases, into an isotropic vortex line liquid. 


Appendix D: Inter-component interaction, and its effect on 
density- and vortex-lattices. 
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FIG. 9: (Color online) Specific heat (top) and helicity moduli 
of both components (bottom), forA^ = 2,/ = 1/16, Cc; = 1.0 
and 7/ = 0.5. Ai ^ « 0.53 we see a clear anomaly in the 
specific heat accompanied by a sharp jump in the 
longitudinal helicity moduli of both components. The 
transverse helicity moduli remain at zero throughout the 
transition. The insets in (a) show the structure factors at the 
high- and low-temperature sides of the transition, 
respectively /3 = 0.528 and P = 0.534. This clearly shows 
that the sharp anomaly in the specific heat separates an 
isotropic phase from a phase with hexagonal order. 


In this appendix, we include more detailed figures of the 
vortex- and density-structures in real- and reciprocal space, as 
the inter-component interaction 2 ( 7 / — cc) |7/i p I 7/2 P is varied, 
to supplement the points made in Sections [IV A| and IV C 


Figs. [m and illustrate how the vortex lattice, and 
the component densities reconstruct as the inter-component 
density-density interaction (7/ — cc)|7/ip|7/2p changes. We 
do this by fixing y at 5.0 and 3.0, respectively, and tuning 
UJ. The inverse temperature, P, is also fixed in both tableaux. 
Common in both figures is that the vortices first form two in¬ 
terlaced square lattices for sufficiently small uj, and, by exten¬ 


sion, large inter-component coupling. Then the lattices recon¬ 
struct into a hexagonal structure. Note that the hexagonal lat¬ 
tices of the two components start out slightly shifted with re¬ 
spect to each other, but becomes completely cocentered when 
UJ > y. This final state corresponds to an attractive inter¬ 
component coupling. 


The behavior of the amplitude densities is explained in 
Sec. |IV A[ and we can compare the reasoning to the top two 
rows of Figs. 11 and 12 First of all, the presence of a vor¬ 
tex locally suppresses the amplitude, which again may affect 
the immediate neighborhood, depending on the value of the 
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inter-component coupling. For strong repulsive couplings a 
suppression of the amplitude of one color in an area leads to 
an enhancement of the amplitude of the other color in the same 
area. The absence of vortices in the neighborhood then leads 
to the opposite effect. This causes staggering of the amplitude 
densities and formation of distinct vortex sublattices. Consid¬ 
ering carefully the range of variation in the amplitudes, it is 
seen that there are rather large gradients for the square struc¬ 
tures. When the coupling is only weakly repulsive or even 
attractive there is a much less dramatic effect. The variations 
in the amplitudes are much smaller; there is little to no stag¬ 
gering. 
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The first column of Figs. and [T^ is in a different class 
from the rest. Here cc = 0, and we are in the SU(2) regime. 
Fig. [^further illustrates the wide variety of ground states ob¬ 
tainable here. This tableau, in contrast to the two previous 
ones, has a fixed 77 and cc, while we vary the inverse temper¬ 
ature /3 from column to column. These pictures are all taken 
from a single simulation, evolved through Monte Carlo sam¬ 
pling from a single randomized initial state as P is increased. 
The vortex lattice initially forms at around /3 = 0.7, and 
evolves continuously. It continues to evolve even at the lowest 
temperatures (P = 1.5) used in the simulation. This pattern is 
common in all simulations done with similar parameter sets. 


FIG. 10: (Color online) Specific heat (top) and helicity 
moduli of both components (bottom), for A^ = 2,/ = 1/16, 
u; = 1.0 and r] = 2.0. At [3 ^ 1.11 we see a clear anomaly in 
the specific heat accompanied by a sharp jump in the 
longitudinal helicity moduli of both components. The 
transverse helicity moduli remain at zero throughout the 
transition, except for the x-directed modulus which drops to 
a negative value at a point well separated from the transition. 
The insets in (a) show the structure factors at the high- and 
low-temperature sides of the transition, respectively 
p = 1.110 and P = 1.112. This clearly shows that the sharp 
anomaly in the specific heat separates an isotropic phase 
from a phase with square order. 
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always arrange themselves in a hexagonal structure, which is 
also seen in the structure factors. 
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FIG. 11: (Color online) Tableau illustrating the different density and vortex lattices in real space, as the parameter uj increases, 
i.e. as the inter component density-density interaction 2 ( 7 ^ — cc) p |V ^2 P decreases. This interaction promotes a square 
density- and vortex lattice. The parameters /, [3, and r] are fixed to / = 1/64, /3 = 1.5 and r] = 5 while uj is increased from 0.0 
to 6.0 horizontally. The six rows show, from top to bottom the amplitude densities of components 1 and 2, the vortex densities 
of components 1 and 2, and the structure factors of components 1 and 2. For cc = 0, which is the SU(2)-symmetric case, the 
system exhibits a dimerized phase in component 1, which is complementary to a dimerized phase in component 2, shifted with 
respect to that of component 1 by an amount corresponding to the lattice constant of the density lattice. The ground state, where 
the roles of components 1 and 2 are inters witched, is degenerate with the illustrated phase. Note that an area of the system with 
a high vortex density always corresponds to an area with a low amplitude density. For the SU(2)-symmetric case, U(l)-vortices 
are not topological. When cc 7 ^ 0, the SU( 2 ) symmetry is broken down to U(l) x U(l), and U(l) vortices are topological. The 
reduction of the interaction 2{r] — uj)\^jji\‘^\'ip 2 \‘^ reduces the tendency towards formation of square density- and vortex-lattices, 
leading to an eventual reconstruction to a standard hexagonal vortex lattice, and hence a hexagonal density lattice. 
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FIG. 12: (Color online) Tableau illustrating the different density and vortex lattices in real space, as the parameter uj increases, 
i.e. as the inter-component density-density interaction 2 ( 7 ^ — Cc;)|'0ip|'02p decreases. This interaction promotes a square 
density- and vortex lattice. The parameters /, [3, and r] are fixed to / = 1/64, /3 = 1.5 and r] = 3 while uj is increased from 0.0 
to 5 horizontally. The six rows show, from top to bottom the amplitude densities of components 1 and 2, the vortex densities of 

components 1 and 2, and the structure factors of components 1 and 2. 
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FIG. 13: (Color online) Tableau illustrating the different SU(2) density and vortex lattices in real space, as the parameter [3 
increases. The parameters /, r], and uj are fixed to / = 1/64, = 1.0 and uj = 0.0 while p is increased from 0.9 to 1.3 

horizontally. The six rows show, from top to bottom the amplitude densities of components 1 and 2, the vortex densities of 
components 1 and 2, and the structure factors of components 1 and 2. Note how the vortex structures and the density structures 
always track, that is, an area of the system with a high vortex density always corresponds to an area with a low amplitude 

density. 
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